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ABSTRACT 

The young stars near the supermassive black hole at the galactic center follow orbits that are 
nearly random in orientation and that have an approximately "thermal" distribution of eccentricities, 
A^(< e) ^ e^. We show that both of these properties are a natural consequence of a few million years' 
interaction with an intermediate-mass black hole (IBH), if the latter's orbit is mildly eccentric and if 
its mass exceeds approximately 1500 solar masses. Producing the most tightly-bound S-stars requires 
an IBH orbit with periastron distance less than about 10 mpc. Our results provide support for a 
model in which the young stars are carried to the galactic center while bound to an IBH, and are 
consistent with the hypothesis that an IBH may still be orbiting within the nuclear star cluster. 
Subject headings: galaxies: active - galaxies: evolution - quasars: general 



1. INTRODUCTION 

Observations of the inner parsec of the Milky Way 
reveal two groups of massive young stars near the su- 
permassive black hole (SBH). A group of roughly 40 
stars move in approximately circular orbits that ex- 
tend inward to a tenth of a parse c from the SBH 
(jPaumard et al. 20061: iLu et al. 20081 ). Another group 
of roughly 20 stars, the S-s tars, follow eccentric or - 
bits well inside the disk stars (lEckart and Genzel 19971 ). 
The orbital periods of the S-stars are as short as 15 
years and the orbits of these stars provide strong con- 
straints on the mass a .nd size of the cen tral dark object 
(jGillessen et al. 20081: [Ghez et al. 20"08l ). 

The presence of young stars so near the galactic center 
is a puzzle, since giant molecular clouds, which are the 
sites of star formation elsewhere in the galaxy, would be 
unabl e to collapse and fragment in the tidal field of the 
SBH H orris 199"3[ ). An orbiting gas cloud would instead 
be sheared into a set of rings; shocks between the rings 
would dissipate energy and the gas would settle into a 
thin disk. If such a disk reaches a critical density, its 
self-gravity can overcome the tidal forces from the SBH 
allowing stars to form (iLevin and Beloborodov 20031 : 
iNavakshin and Cuadra 20051 ). This model has been 
shown to successfully reproduce the observed prop- 
erties of the young disk stars, including their top- 
heavy mass function an d their mildly noncircular orbits 
(|Bonnell and Rich 20081 ). But it fails to account for the 
S-stars, which move on eccentric, randomly-oriented or- 
bits much closer to the SBH. 

A number of models have been proposed to ex- 
plain the S-star s but none is c ompletely satisfactory 
(jAlexander 20051 : iPaumard 2008[ ). The S-stars could 
be old stars that migrated inward and were "rejuve- 
nated" by tidal heating or collisions with other stars. 
However their relativel y normal sp ectra argue against 
such an exotic history (|Figer 2008? ) . Capture of young 
stars onto tightly-bound orbits via three-body exchange 



interactions involving compact remnants (e.g. stellar- 
mass black holes) may e xplain some of the S-stars 
(I Alexander and Livio 20041 ). A similar model assumes 
that the S-stars were originally in binary systems; the 
more massive component of the binary was ejected during 
a close passage to the SBH, leav ing the lower-mass star 
behind (|Gould and Quillen 2003f ). Both of these models 
require an ad hoc reservoir of new stars at large radii, 
as well as some mechanism for placing these stars onto 
plunging orbits soon after their birth, so that they can 
pass near the SBH where the chance of a three-body ex- 
change is appreciable. The models also have difficulty 
reproducing the observed distribution of S-star orbital 
eccentricities. 

An alternative scenario posutlates that the young 
stars formed in a giant molecular cloud, far enough 
from the SBH that tidal forces did no t preclude col- 
lapse and fragmentation (jGerhard 2001h . The newly- 
formed star cluster then migrated inward via dynami- 
cal friction before tidal forces from the SBH dispersed 
it. The presence of an intermediate-mass black hole 
(IBH) at the center of the clu ster would assist in the 
transport (Hansen and Milosavl jevic 200"3l ): in the ab- 
sence of the IBH, the cluster would be completely dis- 
rupted by tidal stresses at a distance of one par - 
sec from the SBH (jPortegies Zwart and McMillan 20M) . 
This model requires a relatively high mass (> IO^Mq) 
and density for the cluster, although not much greater 
than what is observed in existing galactic center 
star cluster s like the Arches and the Quintuplet 
(jFiger 2008f l. While the timescale for formation of 
the IBH is uncertain, simulations suggest that col- 
lisions between stars could form a massive remnant 
in a time shorte r than the time for cluster inspi- 
ral and dissolution (jPortegies Zwart and McMillan 20021 : 
iFreitag et al. 20061 ). This model is also appealing since 
IBHs provide potential solutions to a number of other 
outstanding pr oblems , including the origin of the hyper- 
velocity stars (iLevin 20061 ). the structure of the stel- 
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lar disks (|Leviii et al. 2005t iBerukoff and Hansen 20061) . 
and the growth of SBHs ( Portegies Zwart et al. 2006( 1 

Here we show that the infaUing star cluster model can 
also reproduce the peculiar orbits of the S-stars. 

2. ASSUMPTIONS 

We begin by summarizing two results from re- 
cent TV-body simulat ions of the galactic center 
(iBaumgardt et al. 20061: iMatsubavashi et al. 20081 : 
iLockmann and Baumgardt 2008D . 

• Inspiral of an IBH slows dramatically when it 
reaches a distance ~ 10(g/10~^) mpc from the 
SBH, where q is the ratio of IBH to SBH masses; 
this distance is comparable to the sizes of the S- 
star orbits if g « lO"^, i.e. if Mibh « IO^-^Mq. 
Stalling occurs when the total binding energy in 
background (bulge) stars within the IBH orbit is 
comparable to that of the IBH itself; at this sepa- 
ration, most of the background stars that can ex- 
change energy with the IBH are rapidly removed 
via the gra vitational slingshot an d the frictional 
force drops (jBegelman et al. 1980( ). 

• The orbit of the IBH is likely to be very eccen- 
tric at this late stage, e w 0.5 or greater, the exact 
value depending on the initial orbit of the star clus- 
ter containing the IBH and on the detailed history 
of IBH-star interactions after the cluster has been 
tidally removed. 

We note that all of the A''-body simulations cited above 
assumed an initially steep, p ^ 7--1-4 _ ^-i.75 £jgjj_ 
sity profile around the SBH. In fact, there is evi- 
dence for a "hole" or dip in the density of the domi- 
nant, late- type (o l d) stellar population inside ~ 0.5 p c 
(|Figer et al. 20031 : ISchodel et al. 20071 : IZhu et al. 2008f ). 
Accounting for the observed dip in the stellar densi- 
ties would strengthen both results summarized above: 
a lower density would decrease the dynamical friction 
force on the IBH and lengthen the stalling phase; and a 
shallower density profile is more con ducive to eccentricity 
increases (iGould and Quillen 20031) . 

If sufficiently high eccentricities are reached, e ^ 0.99, 
gravitational radiation losses can shrink the SBH-IBH 
separation from ^ 10 mpc to coalescence in less than 10* 
yr. In what follows, we assume that such extreme eccen- 
tricities are not attained and that the semi-major axis of 
the IBH orbit remains essentially unchanged for times 
comparable to S-star main-sequence lifetimes. How- 
ever we argue below that these assumptions may not 
be strictly necessary for producing the effects that we 
observe. 

Prolonged gravitational interaction with the IBH 
can then scatter the young stars out of the thin 
disk into which they were originally deposited. The 
IBH acts on the stars in much the same way that 
Jupiter acts to scatter comet s in the Solar system 
(|Hansen and Milosavlievic 20031 ) - with the important 
difference that Jupiter's orbit is nearly circular, while 
the IBH orbit is eccentric. 

3. METHODS 

We carried out an extensive set of long-term N- 
body simulations to evaluate the effects of this inter- 
action. Initial conditions consisted of a binary black 
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Fig. 1. — Initial (left) and final (right, after 1 Myr) orbits of stars 
in a simulation with IBH semi-major axis a = 15 mpc, eccentricity 
e = 0.5, IBH-SBH mass ratio q = 0.001, and F = 0.2. Top panels 
show the view looking perpendicular to the IBH orbital plane and 
bottom panels are from a vantage point lying in this plane. The 
IBH orbit is the heavy black curve in all panels; the unit of length 
is milliparsecs. The initially disk-like, co-rotating distribution of 
stars is converted, after 1 Myr, into an approximately isotropic 
distribution of orbits with a range of eccentricities, similar to what 
is observed for the S-stars. 

hole representing the SBH-IBH pair and a set of 19, 
ten-solar-mass stars. Integration of each 21-body sys- 
tem was carried out using the com puter program AR- 
CHAIN (jMikkola and Merritt 20081 ). which incorporates 
an algorithmically regularized chain structure and time- 
transformed l eapfrog to deal with near- collisions between 
the particles (jMikkola fc Merritt 2006f ). Post-Newtonian 
corrections to the equations of motion were included up 
to order PN2.5. Mikkola & Merritt (2008) present a de- 
tailed description of the ARCHAIN algorithm as well as 
results of performance tests when the algorithm is ap- 
plied to galactic center problems very similar to the one 
treated here. 

Initial conditions for the "star" particles were gener- 
ated as follows. 1. The position rapo and velocity Wapo 
of the IBH at apastron were computed. 2. A velocity 
of magnitude -Fwapo was added with random direction 
to the IBH velocity. 3. The Keplerian elements of the 
resulting orbit with respect to the SBH were computed, 
ignoring the presence of the IBH. 4. A random value was 
assigned to the argument of the periastron, and a star 
was placed at a random phase on this orbit. This scheme 
produced an initial distribution of stars about the SBH 
that mimicked the phase-mixed distribution expected for 
a pop ulation of stars that were ti dally stripped from the 
IBH (jBerukoff and Hansen 20061 ) . with a (small) thick- 
ness determined by F (see Fig. 1); all stars were initially 
orbiting in the same sense around the SBH. 

In addition to F, the parameters that were varied were 
the mass ratio q, semi-major axis a, and eccentricity e 
of the SBH-IBH binary. About 300 integrations were 
carried out, each for a time of ~ 5 Myr, based on an 
assumed SBH mass of 4.5 > 
Galactic center of 8.4 kpc. 

4. RESULTS 
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Fig. 2. — Evolution of the Rayleigh parameter R = TZ/N that 
measures the degree of randomness of the stellar orbital orienta- 
tions, in 12 integrations, all with IBH semi-major axis a = 30 mpc. 
Black (solid) lines: q = 10~^; green (dashed) lines: q = 5 X 10~*; 
blue (dotted) lines: q = 2.5 X 10~*; red (dashed) lines: q = 10~*. 
The shaded regions show the 90% confidence bands expected for 
a random, isotropic distribution of orbital orientations of 17 stars 
(the mean number of bound stars in the simulations). Gravita- 
tional perturbations from the IBH produce a nearly random distri- 
bution of orbital orientations after a few Myr, as long as the IBH 
mass and orbital eccentricity are sufficiently large. 




Fig. 3. — Evolution of the distribution of stellar orbital ec- 
centricities e in a set of simulations with IBH orbital parameters 
g = 5 X 10~*, a = 15 mpc, e = 0.5, F = 0.2. Each line is an 
average from six integrations with different random realizations of 
the initial conditions. Six times are shown: t = (thick black 
curve), t = (0.01,0.02,0.04,0.2) Myr (thin blue lines), and t = 1 
Myr (thick red line). The distribution is essentially unchanged at 
times greater than 1 Myr. Dashed line shows a "thermal" distri- 
bution, N cc , and op en circles are the ecce ntricity distribution 
observed for the S-stars IIGillessen et al. 20081 ). 



Figure 1 illustrates the evolution of the stellar orbits 
in one integration with A/jbh = 4500 Af©. In roughly one 
Myr, stars are scattered by the IBH (and occasionally by 
other stars) onto orbits with a wide range of eccentrici- 
ties, semi-major axes and orientations. 

We measured the degree of randomness of the or- 
bital planes u sing the Rayleigh (dipole) statistic TZ 
(| Rayleigh 19 191 ) defined as the length of the resultant 
of the unit vectors U^i = 1....N, where li is perpendic- 
ular to the orbital plane of the ith star. Initially, or- 
bital planes are nearly aligned and TZ ~ N, while for 
a random (isotropic) distribution TZ w N^^^. Figure 2 
shows that TZ reaches values consistent with isotropy in 
a few Myr if the IBH mass is greater than ~ 
and if the eccentricity of the IBH orbit is e « 0.5 or 
greater. Eccentricity of the IBH orbit implies a semi- 
periodic forcing of the stars in a direction perpendicular 
to their initial orbital planes, allowing inclinations to be 
"pumped up" to large values after r epeated encounters 
()I3innev 19811 : lErwin &: Sparke 19991 ). When the orbit of 
the IBH is made nearly circular, e ^ 0.2, stellar orbital 
inclinations were found to remain nearly unchanged over 
these time scales. 

The IBH also induces evolution in the eccentricities 
and energies (semi-major axes) of the stars. Eccentrici- 
ties were found to tend toward a "thermal" distribution 
A^(< e) ~ on a time scale of ~ 0.1 Myr for q > 
2.5 X 10~^ (Figure 3); this is simila r to the distributio n 
that is observed for the S-stars (Gillessen et al. 20081 ). 
In the case of orbital energies, perturbations from near 



encounters with the IBH can either increase or decrease 
them. Up-scattering in energy can continue until a star 
becomes unbound and escapes the system; this typically 
happened for at least one of the stars in each integra- 
tion. Scattering to lower (more tightly bound) energies 
tended to be self-limiting: once a star is transferred to 
an orbit such that its apastron lies inside the perias- 
tron of the IBH, it is "decoupled" dynamically from the 
IBH and its ene rgy tends to remain constant thereafter 
(jWetherill 199 It ). Exceptions were observed only in the 
case F = 0; for such cold initial conditions, the stellar or- 
bits remain highly correlated for long times encouraging 
strong interactions. Producing tightly-bound orbits like 
those of the innermost S-stars, e.g. S2, therefore requires 
an IBH orbit with a periastron distance a(l — e) « 10 
mpc. For IBH orbits satisfying this condition, we were 
able to reproduce approximately the full distribution of 
S-star semi-major axis lengths. However we do not con- 
sider this a crucial test since the observed sample is likely 
to be biased by radius-dependent selection effects. 

We tested the extent to which the randomizing effects 
of the IBH are helped by star-star scattering. We carried 
out additional sets of integrations in which the masses of 
the stars were decreased by factors of ten from IOMq to 
10"'^ Mq. No significant dependence on stellar mass was 
observed, leading us to conclude that interaction with 
the IBH is the dominant mechanism responsible for the 
orbital evolution that we observe. 

5. DISCUSSION 
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We have found that the presence of an IBH orbiting 
within the nuclear star cluster at the center of the Milky 
Way can efficiently randomize the orbits of stars near 
the SBH, converting an initially thin, co-rotating disk 
into a nearly isotropic distribution of stars moving on 
eccentric orbits around the SBH. Randomization of the 
orbital planes requires ~ 1 Myr if the IBH mass exceeds 
~ 1500Mq and if the orbital eccentricity ~ 0.5 or greater. 
Evolution to a "thermal" eccentricity distribution occurs 
on an even shorter time scale. The final distribution 
of stellar semi-major axes depends on the assumed size 
of the IBH orbit, but stars with apastron distances as 
small as the periastron distance of the IBH are naturally 
produced. 

Our simulations contribute only a small piece to the 
bigger unsolved puzzle of the origin of the young stars at 
the galactic center. If models for the genesis of the S-stars 
that invoke an inspiralling IBH are deemed otherwise vi- 
able, our results show how the same models can also nat- 
urally reproduce the random and eccentric character of 
the stellar orbits, and all in a time that is less than stellar 
evolutionary time scales - thus providing an essentially 
complete explanation for the "paradox of youth" of the 
S-stars. 

Although we chose parameters for the IBH such that 
gravitational radiation would not alter its orbit apprecia- 
bly in 5 Myr, the rapidity with which the IBH modifies 
the stellar orbits in our simulations suggests that even 
an IBH on a decaying orbit might be able to randomize 
and "thcrmalize" the S-star distribution before coalesc- 
ing with the SBH. Thus, our model docs not necessarily 
imply that an IBH is present, at the current time, within 
the S-star cluster. But if the IBH is there now, its pres- 
ence might be detected in a number of ways: 

1. The IBH will induc e a motion of the SBH 
(jHansen and Milosavlievic 20 03). Upper limits on the 
astrometric wobble of the radio source Sgr A* are 
so far consistent with the presence of an IBH with 
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mass and semi-major axis in the range considered here 
(jGillessen et al. 2008D . 

2. Stars can remain bound to the IBH if its Hill sphere 
is larger than its tidal disruption sphere; this condition 
is satisfied for SBH-IBH separations greater than ^ 0.05 
mpc. The motion of a star bound to the IBH would 
be the superposition of a Keplerian ellipse around the 
SBH and an additional periodic component due to its 
motion around the IBH; the latter would have a velocity 
amplitude ^ 0.1 — 10 times the IBH orbital velocity and 
an orbital frequency from several hours to a few years, 
potentially accessible to astrometric monitoring. 

3. In favorable circumstances, a near encounter of the 
IBH with a star unbound to it could produce observable 
changes in the star's orbit over month- or year-long time 
scales. 

4. In our simulations with q — 0.001 and e = (0.5, 0.7), 
a fraction 13% of the stars were ejected from the SBH- 
IBH system in 5 Myr. A star ejected at ~ 10'^ km 
s^^ requires ^ 10^ yr to move beyond 0.1 pc imply- 
ing a probability ^ 0.2(A^/10'') of observing an escap- 
ing star at any given time in the Galactic center region, 
where A'^ is the number of stars subject to ejection. One 
S-star in fact appear s to be on an escaping trajectory 
l|Gillessen et al. 2008 f). 

Finally we note that some models for IBH formation 
predict that a large numb er of I BHs mi ght co-exist i n 
the galactic center region (jPortegies Zwart et al. 20061 ). 
If so, the rate of orbital evolution of the young stars might 
be even higher than what we observe in our simulations 
with a single IBH. 
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